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SUMMARY ^ 

Numerical calculation of a three-dimensional turbulent flow of a jet in a 

crossflow using a multiple- time -scale turbulence model is presented. The 
turbulence in the forward region of the jet is In a stronger inequilibrium 
state than that in the wake region of the jet, while the turbulence level 
in the wake region is higher than that in the front region. The calculated 
flow and the concentration fields are in very good agreement with the 
measured data, and it indicates that the turbulent transport of mass, 
concentration and momentum is strongly governed by the Inequilibrium 
turbulence. The capability of the multiple- time -scale turbulence model to 
resolve the inequilibrium turbulence field is also discussed. 
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coefficient for uj -velocity correction 
normalized concentration 

turbulence model constants for €p equation (i-1,3) 
turbulence model constants for e t equation (i-1,3) 
eddy viscosity coefficient 
constant coefficient (-0.09) 
diameter of the circular jet 

index for spatial coordinate (i-1,2,3 and j —1 ,2,3) 
turbulent kinetic energy (k-kp+k t ) 
turbulent kinetic energy in production range 
turbulent kinetic energy in dissipation range 
pressure 

incremental pressure 
production rate 

free-stream velocity of crossflow 

time averaged velocity (-{u,v,w}) 

jet velocity averaged across jet cross-section 

wall coordinate based on friction velocity 

energy transfer rate 

dissipation rate 

molecular viscosity 

molecular diffusivity 

turbulent viscosity 

density 

turbulent Schmidt number 

turbulent Prandtl number for i-equation, i-{kp, ep t k t , 6 t ) 
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INTRODUCTION 


Numerical calculation of a circular jet exhausting Into a cross flow [1] 
is presented. The circular jet in a uniform crossflow is schematically 
shown In Fig. 1. Turbulent flows similar to a jet In a crossflow can be 
found in a number of engineering applications. For example, in gas turbine 
combustors, a number of circumferentially distributed jets are used to 
ensure correct combustion in the flame zone and then to dilute the hot gas 
entering the turbine. Experimental investigations of air jets in crossflows 
have been made to better understand turbulent flows In such engineering 
applications even though air jets in crossflows are by far simpler than 
those in a gas turbine combustor. Compilations of various experimental 
investigations of jets in crossflows can be found in Crabb et al. [1] and 
Khan [ 2 ] . 

With the recent advances in numerical methods to solve the Navier-Stokes 
equations, a number of numerical simulations of jets in crossflows have 
been reported. A compilation of various numerical Investigations of the 
flows can be found in Claus and Vanka [3]. Earlier numerical calculations 
of the flows [4,5] have been devoted to the development and the 
verification of the numerical methods, and only a very small number of grid 
points were used to discretize the entire flow domain due to the limited 
capability of computers. More recently, a realistic number of grid points 
began to be used to calculate the flows using k-e turbulence models [6], 

The numerical results obtained using the relatively fine meshes show 
improved comparison with the measured data in a certain part of the flow 
domain and worse agreement with the measured data in the other part of the 
flow domain. Claus and Vanka [3] carried out a grid independence study of a 
jet in a crossflow using a k-e turbulence model to identify the cause of 
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the deteriorated comparison with the measured data. Claus and Vanka [3] 
showed that the numerical results obtained using a 128x48x48 mesh are not 
significantly different from those obtained using a 256x96x96 mesh and that 
the deteriorated comparison is caused by the inability of the k-e 
turbulence model to describe the complex turbulence field. In each of the 
above numerical simulations, the upstream region of the jet was excluded 
from the Computational domain. However, Andreopoulos [7] showed that the 
circular jet and the crossflow interact strongly with each other at the jet 
exit and that the influence is propagated toward the upstream region of the 
jet. Thus, the deteriorated numerical results can also be caused by the 
numerical models which can not fully account for the strong interaction 
between the jet and the crossflow. In the present numerical study, the 
boundary for the circular jet is located at one diameter upstream of the 
jet exit so that the strong interaction at the jet exit is also accurately 
simulated. 

Numerical results for various complex turbulent flows obtained using 
two -equation turbulence models, algebraic Reynolds stress turbulence models 
(ARSM) and Reynolds stress turbulence models (RSM) show that these 
turbulence models can not accurately describe the turbulence fields of 
various complex turbulent flows [8]. One common inability of the 
two-equation turbulence models, ARSM and RSM is that these turbulence 
models can not account for " inequilibrium turbulence” due to the use of a 
single time scale to describe both the turbulent transport and the 
dissipation of the turbulent kinetic energy. The "inequilibrium turbulence" 
is explained in the "Multiple-Time-scale Turbulence Model” section. The 
multiple- time -scale turbulence model [9] (hereafter abbreviated as the M-S 
turbulence model for convenience) yields accurate numerical results for 
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widely different classes of complex turbulent flows (e.g., turbulent flows 
subjected to extra strains caused by streamline curvatures, interaction of 
multiple number of turbulence fields, and shock wave- turbulent boundary 
layer interactions). The complex turbulent flows to which the present M-S 
turbulence model has been applied as yet include a wall-jet flow, a shear 
layer with wake-boundary layer interaction, a backward- facing step flow, a 
confined coaxial swirling jet, turbulent shear layers over curved surfaces, 
separated transonic turbulent flows over a curved hill and reattaching 
shear layers in a divergent channel. It can be seen in [9-12] that the 
numerical results for these complex turbulent flows obtained using the M-3 
turbulence model are in as good agreement with the measured data as those 
obtained using an optimized k-e, ARSM, or RSM turbulence model for each 
flow case. The capability of the M-S turbulence model to solve widely 
different complex turbulent flows is attributed to its capability to 
resolve the inequilibrium turbulence. This capability is discussed in the 
"Multiple- time-Scale Turbulence Model” section. 

The fluid flow in the near -wall region of the jet exit Is subjected to a 
large mean flow strain rate. The near-wall turbulence field, intensified by 
the large mean flow strain rate, can Influence the entire fluid flow in the 
downstream region of the jet. Thus the near -wall turbulence field In the 
vicinity of the jet exit needs to be resolved accurately In order to 
correctly predict the entire flow field. In the present numerical 
simulation, the near-wall turbulence is described by a "partially low 
Reynolds number" near-wall turbulence model [13], In the model, only the 
turbulent kinetic energy equations are extended to include the near-wall 
low turbulence region and the energy transfer rate and the dissipation rate 
inside the near-wall layer are obtained from algebraic equations. It was 
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shown in numerical calculations of turbulent flows over curved hills [10], 
transonic turbulent flows over a curved surface with shock wave -boundary 
layer interactions [11] , and reattaching shear layers in a divergent 
channel [12] that the "partially low-Reynolds number" near-wall turbulence 
model (when used together with the M-S turbulence model) accurately 
predicts the near-wall turbulence fields. 

The present numerical method is a finite volume method based on a 
pressure correction algorithm [14-16], In the method, all flow variables, 
except pressure and concentration, are located at the same grid points, 
while pressure and concentration are located at the centroid of a cell 
formed by the eight neighboring velocity grid points. The pressure 
correction algorithm is described in the following section. Calculations of 
a three-dimensional lid- driven cavity flow and a laminar flow through a 
90°-bend square duct can be found in [15]. It is shown in [15] that the 
numerical results for the cavity flow obtained using the present numerical 
method compare more favorably with the measured data than those obtained 
using a formally third order accurate quadratic upwind interpolation 
scheme. It Is also shown in [15] that the present method yields a grid 
independent solution for the curved duct flow with a very small number of 
grid points and that the method yields quickly and strongly convergent 
numerical results. Application of the same numerical method for 
two-dimensional flows can be found In [10-12], 


NUMERICAL METHOD 

The incompressible turbulent flow equations are given as; 


d 

(PU4) « 0. (1) 

3x j 
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3 d 3u^ duj 1 3 2 

(pu iU j) (p+M t )( + —)\ (p + — pk) (2) 

3xj [ 3xj dx i J 3x£ 3 

where repeated Indices imply summation over the index unless otherwise 
stated. The convection- diffusion equation for the concentration is given 
as; 

3 3 ' 3 c | 

(pujc) — ) f " 0 (3) 

3xj 3xj a d 3xjJ 

where a d “0.75 is used in the present study. Due to the strong large eddy 
mixing, the molecular diffusivity can be ignored or formally approximated 
as Md'^A'd* and neither of the approximations influence the numerical 
results significantly. 

In the numerical method, the conservation of mass equation is replaced by 
a pressure correction equation given as: 


d (i ap-l a/uj* 

ax j Uj ax jJ ax j 


(4) 


where Uj denotes the velocity which may not satisfy the conservation of 
mass as yet and the last term represents the mass imbalance (see [15] for 
details). The flow equations, the concentration equation, and the 
turbulence equations are solved by a finite volume method using the 
power-law upwind differencing scheme [14]. As all the central-differenced 
finite volume equations for self-adjoint second order elliptic partial 
differential equations are strongly diagonally dominant, the discrete 
pressure correction equation obtained by applying the standard finite 
volume method to eq. (4) is strongly diagonally dominant even for highly 
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skewed and graded meshes. 

For completeness, the velocity-pressure decoupling that occurs when 
various pressure correction algorithms are used for pressure-staggered 
meshes is described briefly below. The use of various pressure correction 
algorithms in pressure -staggered meshes does not yield a diagonally 
dominant system of equations for the incremental pressure. In such a case, 
the mass imbalance at a pressure grid point produces large corrections for 
pressures at adjacent pressure grid points and velocity-pressure decoupling 
occurs [15-16]. For clarity, the differences between the present pressure 
correction algorithm and various other algorithms are summarized below. For 
an orthogonal mesh aligned with cartesian coordinates, the present pressure 
correction algorithm and the "momentum interpolation scheme" of Peric et 
al. [17] and Majumdar [18] yield a 7 -diagonal system of equations for the 
incremental pressure, while various other pressure correction algorithms 
yield a 27 -diagonal system of equations which lacks diagonal dominance 
[15-16]. In the momentum Interpolation scheme, the pressure gradient is 
interpolated differently from the other terms in the discrete momentum 
equation to achieve diagonal dominance. In curvilinear coordinates, the 
present pressure correction algorithm yields a 27-diagonal system of 
equations, while the momentum Interpolation scheme always yields a 
7-diagonal system of equations and it can not account for grid skewness in 
the discrete pressure correction equation in its present form. Also a 
specialized interpolation scheme needs to be adopted to obtain a unique 
solution that does not depend on the under -relaxation parameter [18]. On 
the other hand, the present pressure correction method yields a unique 
solution since the incremental pressure is driven only by the mass 
imbalance as shown In eq. (4). 
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MULTIPLE -TIME -SCALE TURBULENCE MODEL 


The anisotropy of the turbulence is the most easily detectable phenomenon 
in a measurement of a turbulent flow. Thus, it was conceived that the poor 
capability of the two- equation turbulence models to resolve complex 
turbulent flows is attributed to the inability of the turbulence models to 
take into account of the anisotropy of the turbulence. Thus the emphasis is 
laid upon improving the ARSM and the RSM. However, a number of numerical 
investigations carried out during the last one and half decades show that 
the ARSM and RSM still can not accurately predict the turbulence phenomena 
occurring in various flows unless the pressure-strain rate correlation is 
optimized for each flow [11-12]. 

A careful examination of semi -empirical data (theoretically derived data 
from a set of measured data) reveals that the inequilibrium turbulence also 
dictates the developments of the mean flow field and the turbulence field. 
Here, the "inequilibrium turbulence” represents the state of a turbulence 
field in which Pr/e t varies rapidly in space so that the shape and the 
frequency domain of the spectral density varies widely in space. The 
spectral density curves shown in Fig. 2- (a) are constructed based on the 
measured data of Klebanoff [19] and Wygnanski and Fiedler [20], It can also 
be seen in Fig. 2- (a) that the generation of the energy containing large 
eddies by the instability of the mean fluid motion occurs in the low 
frequency region and that the peak of the spectral density moves toward the 
high frequency region as Pr/e t is decreased. The semi-empirical c^ for a 
plane jet obtained by Rodi [21] is shown in Fig. 3. It can be seen In the 
figure that c^ is decreased as Pr/e t is increased, and c^ is increased as 
Pr/e t Is decreased. Thus, the developments of the mean fluid flow and the 
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turbulence field are influenced by the spatially varying turbulent eddy 
viscosity, and the spatially varying turbulent eddy viscosity depends on 
the level of the local inequi librium turbulence. Consider k-e turbulence 
models for which the eddy viscosity is given as . Due to the 

use of a constant c^f, the eddy viscosity in a strongly turbulent region 
(Pr/e t >l) is over-predicted and that in weakly turbulent region (Pr~0) is 
under-predicted. As an example, the under-predicted reattachment location 
of a backward- facing step flow obtained using a k-e turbulence model is 
caused by the over-predicted turbulent viscosity along the reattaching 
shear layer [12], 

The variation of c^ as a function of Pr/e t was incorporated into k-e 
turbulence models in the form of "generalized algebraic stress turbulence 
models." The c^ curves by Launder [22] and Kim and Chen [23] are shown in 
Fig. 3. The generalized algebraic stress turbulence models yield accurate 
numerical results for shear layers when used in boundary layer flow 
solvers. However, the use of these turbulence models in elliptic 
(two-dimensional) flow solvers does not easily yield a converged solution 
due to a severe interpolation used in the c^ function. Furthermore, it is 
not clear if the generalized algebraic stress turbulence models can resolve 
the inequilibrium turbulence as cleanly as the M-S turbulence models can, 
since the generalized algebraic stress turbulence models lack many features 
of the M-S turbulence models to be described below. 

The M-S turbulence models [9,24] appeared as a consequence of recognizing 
the inability of various single-time-scale turbulence models (k-e, ARSM, 
and RSM) to accurately describe complex turbulent flows. The 
convection- diffusion equations of the M-S turbulence models were 
established by Hanjelic et al. [24]. The convection-diffusion equations 
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most naturally describe the physically observed turbulence phenomena in the 
sense that the turbulent transport of mass and momentum is described using 
the time -scale of large eddies and the dissipation rate is described using 
the time-scale of fine-scale eddies. Later, Kim and Chen [9] established 
the general form of the load functions based on a physical dimensional 
analysis. The differences between the present M-S turbulence model and that 
of Hanjelic et al . [24] can be found in [11-12], and hence, these are not 
repeated here. The capability of the M-S turbulence model to solve complex 
turbulent flows was further enhanced by incorporating a "partially 
low-Reynolds number" near-wall turbulence model into the M-S turbulence 
model [10-13]. Calculations of widely different classes of complex 
turbulent flows showed that the M-S turbulence model can resolve the 
inequilibrium turbulence and can model the cascade of turbulent kinetic 
energy. These capabilities of the M-S turbulence model are described below. 

The M-S turbulence equations are given below for completeness. The 
turbulent kinetic energy and the energy transfer rate equations for energy 
containing large eddies are given as; 
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where the production rate is given as; 
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The turbulent kinetic energy and the dissipation rate equations for fine 
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scale eddies are given as: 
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and the eddy viscosity is given as; 


k2 

Pt ~ P c fif— (9) 

e P 

The turbulence model constants are given as; oj^-0.75, < 7 j ct - 0 . 75 , a €p -1.15, 
a ct~ 1-15, Cp^— 0.21, Cp 2 ~l . 24 , c p 3 “ 1.84, c^2_”0-29, c^ 2 ” 1.28, and c^ 2 “l- 86 . 

The capability of the M-S turbulence model to resolve the inequilibrium 
turbulence depends largely on the load functions of the e p and e t equations 
and the way the turbulence model constants are established. The load 
functions of the e p and £ t equations are obtained from a physical 
dimensional analysis [9], and the establishment of the model constants are 
based on the assumptions that the turbulence field of a uniformly sheared 
flow can approach an asymptotic state in which Pr/£ t becomes a constant and 
that the ratio of f^Ap depends on the ratio of Pr/£ t . The first assumption 
that such an asymptotic state can exist was shown by Harris et al. [25], 
and later, was confirmed by Tavoularis and Karnik [26], In such asymptotic 
states of uniformly sheared flows, the diffusion term vanishes, and the 
asymptotic ratio of kp/k t can be obtained by dividing eq. ( 5 ) by eq. ( 7 ), 
i.e. , 


kp Dkp Pr/e t - tp/c^ 

k t Dk t £ p/ £ t * 1 


( 10 ) 
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It can be seen in eq. (10) that the existence of the asymptotic ratio of 
kp/k t depends on the realizability of the second assumption. The second 
assumption that the ratio of € t Ap depends on the ratio of Pr/e t can be 
verified by numerical results posterily, or it can be verified indirectly 
by comparing the M-S eddy viscosity equation with that of the generalized 
algebraic stress turbulence models shown in Fig. 3. The eddy viscosity, eq. 
(9) , can be rewritten in a form compatible with that of the generalized 
algebraic stress turbulence models, I.e., 

k 2 

l*t ~ (11) 

e t 

where c^-c^f (c t Ap) is the eddy viscosity coefficient, the variation of c^ 
is described by the ratio of £ t Ap» and e t Ap I s a function of P r A t . Thus 
the second assumption can be justified within the context of the 
generalized algebraic stress turbulence models. 

The three inequilibrium turbulence levels (A, B, and C) imbedded into the 
M-S turbulence model (or used in determining the turbulence model 
constants) are also shown in Fig. 3. The measured data that corresponds to 
the point A in Fig. 3 (i.e., Pr/c t -1.5) can be found in Harris et al. [25] 
and in Tavoularis and Karnik [26]. The value of € t A p -0.95 can be estimated 
from Fig. 3. The ratio of kp/k t for the data point A is obtained to be 9.0 
from eq. (10). For turbulent flows in an equilibrium state (point B in Fig. 
3), Pr-e t , and e p has to be equal to both of them to maintain the 
equilibrium state. In this case, eq . (10) becomes indeterminate; and the 
ratio of kp/k t ==4.0 can be obtained from a near-wall analysis of turbulent 
flows In equilibrium state [9]. In the free stream region of turbulent 
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flows, the production rate vanishes. Such a case is represented by Point C 
in Fig. 3 and the ratio of e t /cp-2.5 can be estimated from the generalized 
algebraic stress turbulence model of Kim and Chen [23], The ratio of 
kp/k t ~0.7 can be obtained from eq (10). 

The three ratios of kp/k t obtained in the above analysis show that 


V 

> 

V 

> 

V 


A 


B 

- k t. 


where A, B and C denote each turbulence state marked in Fig. 3. The 
implication of the above inequality is illustrated in Fig. 2-(b), where the 
frequency domain is divided into two parts by the simplified split-spectrum 
[9,24], The spectral density A in Fig. 2-(b) belongs to a production 
dominated region (P r /£ t >l) , B belongs to a near equilibrium region 
(P r /e t ssl), and C belongs to a region where production rate vanishes (i.e., 
free stream region). The spectral density curves A, B, and C also describe 
the cascade process of turbulent kinetic energy that the energy containing 
large eddies (characterized by low frequencies) are generated by the 
instability of the mean fluid flow, the large eddies become finer eddies 
(characterized by higher frequencies), and the fine scale eddies are 
dissipated by the viscous force. It can be seen in Fig. 2-(b) that kp/k t is 
greater for larger eddies and kp/k^ becomes smaller for finer-scale eddies. 
Thus the ratio of kp/k t is determined by the shape and the frequency domain 
of each spectral density. 

Calculations of various shear layers (boundary layer flows) using the M-S 
turbulence model always reproduce the Imbedded inequilibrium turbulence 
states. For highly complex turbulent flows, large eddies generated in the 
upstream region are convected in the downstream direction. In this 
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downstream region, the relationship between P r /e t and c t / c p influenced 
by the convected eddies (i.e., the large value of kp/k t ) and the numerical 
results exhibit the trend of imbedded inequilibrium turbulence. 

NUMERICAL RESULTS 

The flow domain of the jet in a uniform crossflow [1] is shown in Fig. 1, 
The jet velocity averaged across the cross-sectional area of the pipe is 
27.6 m/sec, the free stream velocity of the cross flow is 12 m /sec, and the 
diameter of the jet is 0.0254 ra. In the experiment [1], the concentration 
field was measured by injecting helium gas (He) into the circular jet. The 
concentration of the helium is one percent of the air-helium mixture at the 
jet exit, and hence the concentration equation, eq. (3), is solved 
uncoupled from the momentum equations. The symmetric half of the flow 
domain Is discretized by 148x61x94 grid points in x- , y- , and z-coordinate 
directions, respectively. The body- fitted grid near the jet exit is shown 
in Fig. 4. The smallest mesh size in the direction normal to the wall is 
0.6x10“^ m (y + =1.5 based on the fully developed pipe flow) and this mesh 
size is sufficiently small to resolve the near-wall turbulence field in the 
vicinity of the jet exit. The largest mesh size used near the far field 
boundaries is approximately half of the jet diameter. 

The inlet boundary conditions for the tangential velocity, the turbulent 
kinetic energies, and the dissipation rates (€p and c t ) are obtained from 
measured data for a fully developed boundary layer flow over a flat plate 
[19]. The non-dimensional velocity and the turbulent kinetic energy 
profiles are scaled to yield a boundary layer thickness of 0.005 m at the 
inlet boundary. The no- slip boundary condition for velocities, vanishing 
gradient for concentration, and vanishing turbulent kinetic energy are 
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prescribed at the solid wall boundary. A vanishing gradient boundary 
condition is used for all flow variables at the exit boundary. The 
symmetric boundary condition is used on the south (y-0) boundary, and the 
free stream boundary condition is used on the north (y-3.1D) and the top 
(z-7 . 5D) boundaries . 

The bottom boundary is located at one jet diameter upstream of the jet 
exit (z«-lD), and a fully developed pipe flow and a constant concentration 
(c-1.0) conditions are prescribed as the jet inlet boundary. These boundary 
conditions can greatly reduce the uncertainty that can be caused by 
ill -posed boundary conditions at the jet inlet as discussed in the 
following paragraph. The use of the free stream boundary conditions on the 
north and the top boundaries is not a good approximation of the actual 
fluid flow unless these boundaries are located sufficiently far away from 
the jet exit. However, it can be seen in the following that the numerical 
results near the jet exit are not influenced too much by the far field 
boundary conditions [10], The partition between the near-wall layer and the 
external region is located at y + ^100 (based on the fully developed flat 
plate flow) and 10 grid points are allocated inside the near-wall layer. 

The converged solutions are obtained in approximately 1200 iterations, and 
the relative mass and concentration imbalances are 2.5x10“^ and 9.5x10'^, 
respectively . 

The contour plots of the jet velocity, the pressure, and the total 
pressure at the jet exit are shown in Fig. 5, where the increments between 
the contour lines are the same for each contour plot. It can be seen in the 
figure that the jet velocity, the static pressure, and the total pressure 
vary widely across the cross-section. It does not seem possible to 
prescribe a correct boundary condition for the jet if the bottom boundary 
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is located at the jet exit. For example, in previous numerical calculations 
of jets in crossflows, either a constant vertical velocity or a constant 
total pressure was prescribed at the jet exit [2,4], However, the present 
numerical results show that a significant amount of uncertainty can be 
caused by the use of either of these boundary conditions at the jet exit. 

The calculated velocity vectors, pressure, turbulent kinetic energy, 
Pr/e t , € t/cp anc * kpAt are shown in Fig. 6. The velocity vector and the 
pressure contour plot show that the crossflow is decelerated rapidly by the 
jet and thus the pressure is increased in the forward region of the jet. 
Otherwise, the velocity vector and the pressure coritour plot do not show 
that any significant phenomena occur in the forward region. However, the 

complex Pr/e t and contours show that the turbulence field is 

experiencing an enormous evolution in the forward region. These contour 

plots show that the inequilibrium turbulence becomes stronger as a fluid 
particle approaches the jet and that the peak inequilibrium state occurs 
along the interface of the jet and the crossflow. The turbulent kinetic 
energy, the production rate, the enefgy transfer rate and the dissipation 
rate in the wake region of the jet are by far greater than those in the 
forward region of the jet. However, the turbulence in the forward region is 
in a stronger inequilibrium state than that in the wake region. These 
results indicate that the strength of inequilibriura turbulence does not 
necessarily depend on the turbulence intensity. It can be seen from eq. 

(11) that the rapidly varying e t Ap in the forward region of the jet will 
influence the fluid flow significantly. It takes a while for large eddies 
to cascade to smaller eddies. The large ratio of kp/kt * n t * ie wake region 
of the jet Is caused by the large eddies convected from the upstream region 
and those generated in the wake region, see Fig. 6-(a). 
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The calculated vertical velocity profiles in the vicinity of the jet exit 
are compared with the measured data in Fig. 7. The slightly higher 
w-velocity for x/D> 2 is caused by the north boundary which is not located 
far enough from the jet exit. Due to the prescribed velocity profile (which 
is the same as that of the inlet plane) on the north boundary, the jet can 
not expand freely into the north direction and thus the excess mass flow 
rate along the symmetry plane cause the over-predicted w-velocity in the 
region. 

It Is shown in Fig. 8 that the calculated tangential velocity profiles 
along the x-coordinate direction on the symmetry plane are In good 
agreement with the measured data. The reversed flow behind the jet indicate 
that the development of the tangential velocity along the crossflow 
direction is similar to that of the flow over a circular cylinder. However, 
the u- velocity in front of the jet is not brought to zero due to the 
compliance of the jet. 

The u-velocity profiles at four downstream locations are shown in Fig. 9. 
The complexity of the u-velocity profiles are caused by the deflected jet 
and the separated crossflow behind the jet. It can be seen In the figure 
that the calculated results are In good agreement with the measured data 
qualitatively and quantitatively. 

The calculated turbulent kinetic energy distribution along the x-axis of 
the symmetry plane at z/D-0.75 Is compared with the measured data in Fig. 

10. It can be seen In the figure that the trend of the turbulent kinetic 
energy distribution is in excellent agreement with the measured data even 
though the turbulence intensity is under -predicted. 

The calculated concentration profiles at three downstream locations on 
the symmetry plane are shown In Fig. 11, and the concentration contour plot 
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at x/d-8 is shown in Fig. 12- (b). The shape and the peak locations of the 
calculated concentration profiles are in very good agreement with the 
measured data. The slightly smaller magnitude of the concentration is 
caused by the coarse grid inaccuracy in the far downstream region. For 
example, the concentration level near (x/D, y/D, z/D)~(8.0, 0.7, 2.1) shown 
in Fig. 12- (b) is slightly higher than the measured data, and hence the 
concentration level on the symmetry plane is slightly under-predicted. 

The calculated u- velocity, concentration, and turbulent kinetic energy 
contours are shown in Figs. 12-(a), (b) , and (c) , respectively. It can be 
found in Crabb et al. [1] that the present numerical results are in 
excellent agreement with the experimentally obtained contour plots. The 
slight difference between the calculated and the measured u-velocity 
contour plots in the vicinity of z/D»*4.5 is again attributed to the coarse 
grid inaccuracy in the region. Note that the peak concentration occurs in 
the region where u-velocity is minimum. This trend indicate that the 
turbulent transport of the concentration is significantly different from 
that of the momentum. The concentration contour plot exhibits a strong 
similarity with the turbulent kinetic energy contour plot. This trend is 
also In excellent agreement with the experimentally observed distributions 
of the concentration and the turbulent kinetic energy [1] . The excellent 
agreement between the calculated and the measured contour plots indicates 
that the M-S turbulence model can correctly resolve the turbulent 
transports of mass and momentum and that the turbulent transport of mass 
and momentum depends strongly on the inequilibrium turbulence. It is not 
clear as yet if any single- time-scale turbulence model can correctly 
resolve the concentration field as the H-S turbulence model can. 

The three-dimensional particle trajectories are shown in Fig. 13. It can 
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be seen in the figure that the fluid particles passing near the jet exit 
are most easily entrained to the jet. It is also shown in the figure that 
the fluid particles near the jet edge carry less momentum and hence these 
particles are quickly entrained to the helical vortices in the wake region 
of the jet. The particle trajectories show that the large eddy mixing 
occurs in the wide region of the jet edge and that the fluid particles in 
the center region of the jet does not mix easily with the crossflow. The 
concentration profiles shown in Fig. 12 also indicate the same trend of the 
large eddy mixing. 


CONCLUSIONS AND DISCUSSION 

It has been shown that a strong inequilibrium turbulence field is 
characterized by the shape and the frequency domain of the spectral density 
that varies widely in space. The influence of the Inequilibrium turbulence 
on the development of the mean fluid flow (and consequently, on the 
development of the turbulence field itself) can be sensed only through 
semi -empirical data. Thus the influence of the inequilibrium turbulence on 
the mean fluid flow is more difficult to recognize than other turbulence 
phenomena such as the anisotropy of turbulence. The semi -empirical data 
show that the eddy viscosity coefficient becomes smaller in the production 
dominated region, and becomes larger In the dissipation dominated region. 

In the multiple- time -scale turbulence model, the dependence of the eddy 
viscosity coefficient on the Inequilibrium turbulence is reflected In the 
ratio of the dissipation rate (e t ) to the energy transfer rate (€p). In the 
simplified split- spectrum case, the measured spectral density curves show 
that the ratio of kp/k^ Is greater for larger eddies, and becomes smaller 
as the large eddies cascade down to smaller eddies. The calculated kp/k*. 
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using the multiple- time -scale turbulence model also shows the same behavior 
as that observed in experiments. The accurate numerical results for a wide 
class of complex turbulent flows (including the present jet in the 
crossflow) obtained using the present multiple- time-scale turbulence model 
indicate that the turbulent transport of mass, concentration, and momentum 
depends strongly on the inequilibrium turbulence and that the 
multiple- time-scale turbulence model correctly resolves the inequilibrium 
turbulence phenomena. It Is not clear if single- time -scale turbulence 
models can resolve such inequilibrium turbulence phenomena as yet. 

Numerical results for the three-dimensional turbulent flow of a circular 
jet In a crossflow show that the jet and the crossflow interact very 
strongly with each other In the forward region of the jet and that the 
interaction creates a strong inequilibrium turbulence field in the forward 
region of the jet. The strong Interaction between the jet and the crossflow 
at the jet exit also influences the fluid flow and the turbulence field in 
the upstream region of the jet. This results suggests that the upstream 
region of the circular jet needs to be included into the computational 
domain in order to obtain accurate numerical results or to assess the 
predictive capability of a turbulence model. 

The calculated velocity, concentration, and turbulence fields are in good 
agreement with the measured data. Both the calculated results and the 
measured data show that the jet in crossflow is characterized by highly 
complex velocity, concentration, and turbulence fields that are not usually 
found In many other turbulent flows. It is discussed in Crabb et al. [1] 
that the weak vortex shedding does not influence the mean fluid flow 
significantly. The good comparison between the numerical results and the 
measured data is also in agreement with such an observation. The calculated 
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tangential velocity, concentration, and turbulent kinetic energy contours 
at a downstream location show that the peak concentration occurs where the 
tangential velocity becomes local minimum and that the concentration field 
exhibits a close resemblance to the turbulence field. These contour plots 
are in excellent agreement with the measured contour plots. 
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7.5 D 



T 


Figure 1 — Nomenclature and computational domain tor a 
circular jet in crossflow. 



log^Kjm- 1 ] 


(a) Spectral density for inequilibrium turbulent flows, 

A: maximum shear location in a circular jet [20], 

B: center of a circular jet [20], C: free stream region 
of a boundary layer flow in zero pressure gradient [19]. 



K 

k P-to EdK - k ><-K 1 EdK 

(b) kp/k t for Inequilibrium turbulent flow, A: Pr/e t > 1, 
B: Pr/€ t « 1, C: Pr*0.0. 


Figure 2— Spectral density. 
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u/TJ, 




x/D 


x/D 


(a) z/D = 075. 


(b) z/D = 1.35. 


Figure 8.— Tangential velocity (u) profiles on the symmetry plane. 



Figure 9— Tangential velocity (u) profiles 
on the symmetry plane, (a) x/D = 0, 

(b) x/D = 0.5, (c) x/D - 0.75, (d) x/D - 
1 . 0 . 



x/D 

Figure 10.— Turbulent Kinetic energy 
profile on the symmetry plane at 
z/D = 0.75. 



Figure 1 1. — Concentration profiles at 
downstream locations on the symmetry 
plane, (a) x/D = 4, (b) x/D = 6, (c) x/D = 
8 . 
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(a) u-velocity. (b) Concentration, (c) Turbulent 

kinetic 

energy. 

Figure 12. — Contour plots of u-velocity, concentration, 
and turbulent kinetic energy. 



Figure 13. — Three-dimensional particle trajectories of a jet in a crossflow. 
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